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ABSTRACT 

We model the interaction between the wind from a newly formed rapidly rotat- 
ing magnetar and the surrounding supernova shock and host star. The dynamics is 
modeled using the two-dimensional, axisymmetric thin-shell equations. In the first 
~ 10 — 100 seconds after core collapse the magnetar inflates a bubble of plasma and 
magnetic fields behind the supernova shock. The bubble expands asymmetrically be- 
cause of the pinching effect of the toroidal magnetic field, even if the host star is 
spherically symmetric, just as in the analogous problem of the evolution of pulsar 
wind nebulae. The degree of asymmetry depends on E mag / E totl the ratio of the mag- 
netic energy to the total energy in the bubble. The correct value of E mag /E to t is 
uncertain because of uncertainties in the conversion of magnetic energy into kinetic 
energy at large radii in relativistic winds; we argue, however, that bubbles inflated 
by newly formed magnetars are likely to be significantly more magnetized than their 
pulsar counterparts. We show that for a ratio of magnetic to total power supplied by 
the central magnetar E mag /E to t <J 0.1 the bubble expands relatively spherically. For 

E m ag/Etot J> 0.3, however, most of the pressure in the bubble is exerted close to the ro- 
tation axis, driving a collimated outflow out through the host star. This can account for 
the collimation inferred from observations of long-duration gamma-ray bursts (GRBs). 
Outflows from magnetars become increasingly magnetically dominated at late times, 
due to the decrease in neutrino-driven mass loss as the young neutron star cools. We 
thus suggest that the magnetar-driven bubble initially expands relatively spherically, 
enhancing the energy of the associated supernova, while at late times it becomes pro- 
gressively more collimated, producing the GRB. The same physical processes may 
operate in more modestly rotating neutron stars to produce asymmetric supernovae 
and lower energy transients such as X-ray flashes. 

Key words: Stars: neutron; stars: supernovae: general; gamma-rays: bursts; stars: 
winds, outflows; magnetic field; MHD 



1 INTRODUCTION 

In the first few seconds after core collapse in a massive 
star, a proto neutron star (PNS) cools and contracts on its 
Kelvin-Helmholz cooling timescale (~ 10 — 100 s), radiating 
its gravitational binding energy (~ 10 53 ergs) in neutrinos 
l|Burrows!ll986l , IPons et alj|l999h . The cooling epoch is ac- 
companied by mass loss driven by neutrino energy deposi- 
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tion in the atmosphere of the PNS (e.g. jDuncan et al.l 1 19861 . 
Qian fc Wooslevlll99ej . iThompson et al.ll200ll ). 



A subset of neutron stars are inferred to have mag- 

net ic field strengths of ~ 10 14 — 10 15 G ("magnetars"; 

see IWoods fc Thompson! |2004| for a review) . If some mag- 
netars are born with milli s econd rotation period s (e.g ., 
iDuncan fc Thompson! Il992l . IThompson fc Duncan! Il99ah 
the combination of rapid rotation and strong magnetic fields 
makes the winds from young magnetars significantly more 
energetic than the thermally driven winds from slowly rotat- 
ing neutron stars. In addition, as the neutrino-driven mass 
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loss decreases during the Kelvin-Helmholtz cooling epoch, 
the wind from a proto-magnetar becomes increasingly 
magnetically-dominated and the flow eventually becomes 
relativistic. For this reason, proto-magnetars have been con- 
sidered as a possible central engine for long-duration gamma 
ray b ursts (GRBsfl a nd hyper-energetic supernovae llUsavl 



19921, lThompsonlll994 IWheeler et allboool . iThompson et al 



2004 ) , and as a possible source of ultra-high energy cosmic 



rays (jBlasi et al.ll2000l . lAronsH2005t ). 

The discovery that GRBs are at cosmological distances 
confirmed that the isotropic energy scale for the gamma- 
ray emission from GRB s is ~ 10 52 — 10 53 ergs (see, e.g 



IWooslev fc Blo om 2006). However, the interpretation of af- 
terglow observations ("jet breaks") suggested that GRBs 
are powered by collimated jets and that the i ntrinsic en- 
ergy in relativistic material is ~ 10 51 ergs (e.g.. iFrail et al.l 
l200ll ). This interpretation has become less clear in recent 
years because of the complex time-dependence in SWIFT 
X-ray afterglow observations and the lack of evidence for X- 
ray jet breaks in the first ~ 10 days fe.g. JSato et alJl2007L 
iBurrows fc Racusinl l2007h . Nonetheless, the case for colli- 
mated outflows from GRBs is compelling. Theoretically, the 
association of many long -duration GRBs with supernovae 
l|Wooslev fc Blo om 2006) sets the natural energy scale for 
GRBs at ~ 10 jl — 10 52 ergs. In addition, estimates of the 
energy in relativistic outflows in GRBs from late time radio 
observations provide lower limits of the same order, although 
the true energy could in principle be much higher (see, e.g., 
lEichler fc Waxmanll2005h . 

In the collapsar model (e.g., iMacFadven fc Wooslevl 

1999), the collimated outflows from GRBs are accounted 
for by jets produced by an accretion flow onto a cen- 
tral black hole. In the magnetar model, the origin of 
such collimated outflows is less clear. Relativistic magne- 
tized outflows by themselves do not efficiently self-collimate 
(e.g., iLvubarskv fc Eichlerj|200lh . Although observations of 
pulsar wind nebulae (PWNe)- which are intrinsically far 
more relativistic than GRBs - show iet-like features (e.g., 
Weisskopf et all I2000L iPavlov et all |200ll iGaensler etafl 



2002), these are believed to be only mildly relativistic out 



flows produced by the interaction between the pulsar wind 
and the surrounding expanding supernova (SN) remnant 



|Komissarov fc Lvubajskvll2004 iDel Zanna et al.ll2004T ). In 
this paper, we explore the hypothesis that collimated out- 
flows from newly formed magnetars can likewise be produced 
by the interaction between the magnetar wind and the sur- 
rounding host star. 

Our physical picture is that the fast trans-magnetosonic 
magnetar wind shocks on the relatively slow outgoing SN 
envelope, creating a subsonic bubble of plasma and mag- 
netic fields inside its host star. Because of the strong toroidal 
magnetic field and the accompanying pinch, an anisotropic 
pressure distribution between the pole and equator is set 
up within the cavity defined by the SN shock and the in- 
coming magnetar wind. For simplicity we assume that (1) 
an outgoing SN shock has created a central evacuated cav- 
ity and (2) the surrounding host star is spherically sym- 
metric. Assumption (1) allows us to model the problem of 
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interest as a free magnetar wind interacting with the ex- 
panding envelope created by a SN shock that is in turn 
sweeping through the host star. Spectral modeling of the 
hyper-energetic supernovae associated with several GRBs 
suggests massive pr ogenitor stars (e.g.. llwamoto et aflll998l ; 
iMazzali et al.ll2006h . This has been interpreted as indicating 
that GRBs are associated with the formation of black holes. 
However, there is increasing evidence that some Galactic 
magnetars ari se from massive s tars with ZAMS masses of 
&40M Q fe.g. jMuno et alj|200d ). Thus our assumption of a 
successful core-collapse SN leaving behind a rapidly rotat- 
ing magnetar is quite reasonable given current observational 
constraints on the progenitors of magnetars and GRBs. Our 
assumption (2) that the host star is spherically symmetric 
may be conservative. Multi-dimensional simulations of core- 
collapse in the presence of rapid rotation and strong poloidal 
magnetic fields find that the explo sion may occur preferen- 
tially along the rota tion axis (e.g., iLeBlanc fc Wilsonlll970l . 
IBurrows et al J 12007). It is presumably easier to produce a 
late-time collimated outflow in this case, since a low pres- 
sure, low inertia channel has already been created. 

A full magnetohydrodynamic (MHD) simulation of the 
interaction between a magnetar wind and its host star would 
require resolving a very wide radial dynamic range. In addi- 
tion, the physical conditions in the wind at large distances 
in particular, the magnetization of the wind - are not 
fully understood (§2). For these reasons, we believe that it 
is fruitful to solve a model problem that allows one to readily 
explore the parameter space of magnetar-host star interac- 
tions - the thin-shell approximation provides such a model. 

In the thin-shell approximation, one assumes that the 
material swept-up by the wind from the central object is 
confined to a geom etrically thin s hell, whose dynamics is 
then evolved (e.g., iGiulianil |l982| ). This model has been 
extensively applied in the study of the interaction of stel- 
lar winds with their surrounding environment , both in the 
case of momentum driven wi nds (see, e.g., ICantol Il980l ; 
ICanto et al. I ll996l ; IWilkinlk oOO) and i n the case of pressure 
driven winds (e.g., Chevalier fc Luol IT994') . The evolution 
of ma gnetized PWNe boun ded by an expanding SN rem- 
nant (Bcgclma n fc Lilll992l ) is the closest analogue to the 
problem we consider in this paper. In a number of cases, 
more detailed numerical simulations ha ve confirmed the va- 
lidity of the thin- shell model (see, e.g.. IStevens et al]|l992l ; 
iBucciantinil |2003 for hydrodynamical examples). Most im- 
portantly for our purposes, axi s ymm etric relativistic MHD 
simula tions by Ivan der Swaluwl (20031 ) andlPef Zanna et al.l 
(2004) have shown that the overall shape o f PWNe resem- 
bles th at predicted by the thin-shell model of Bcgclman & Li 
( 1992). For these reasons we believe that the thin-shell shell 
approximation is a useful tool for studying the structure and 
evolution of bubbles formed by magnetar winds inside their 
progenitor stars. In addition, these calculations can define 
the most interesting parameter space for future relativistic 
MHD simulations. 

The remainder of this paper is organized as follows. 
In § 2 we discuss the general properties of proto-magnetar 
winds, and how they evolve in the ~ 100 seconds after core- 
collapse. We also discuss the equilibrium structure of the 
magnetized bubble created by the magnetar wind behind the 
SN shock. Section 3 summarizes the thin-shell equations. In 
§ 4 we present our results for the evolution of the SN shock 
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Figure 1. Evolution of a magnetar wind as a function of time 
since core bo u nce, based on the evolutionary calculations of 
iMetzeer et al.l {2007). Top: Energy loss rate Bottom: Magneti- 
zation at the light cylinder. Models are for an initial period of 
P = 1 ms, B = 10 15 G (dot dashed), 6 = 3x 10 15 G (solid), 
and B = 10 16 G (dotted) and an initial period of P = 2 ms and 
B = 10 16 G (dashed). 



due to the asymmetric pressure produced by the interior 
magnetized bubble. In § 5 we summarize our conclusions 
and discuss the implications of our results for understand- 
ing observations of long-duration gamma-ray bursts, X-ray 
flashes, and asymmetric supernovae. In the Appendix we 
present self-similar solutions that provide insight into how 
the shape of the bubble is related to its magnetization and 
the conditions in the ambient medium. 



2 PROTOMAGNETAR EVOLUTION AND 
BUBBLE STRUCTURE 

2.1 Protomagnetar Wind Evolution 

In order to set the stage for the thin-shell evolutionary calcu- 
lations that follow, this section summarizes some of the key 
prop erties of outflows from young rapidly rotating magne- 
tars. IMetzeer et alJ (|2007h (hereafter M07) solved the one- 
dimensional (monopolar) neutrino-heated non-relativistic 



MHD wind problem for conditions appropriate to young 
magnetars. These calculations provide the mass-loss rate 
(M) and energy-loss rate (E) from the magnetar as a func- 
tion of parameters such as the neutrino luminosity, mag- 
netic field strength B, and rotation rate Q. The calculation 
of M is applicable even if the wind is relativistic because the 
mass-loss is set in the non-relativistic subsonic portion of the 
wind at small radii. The calculations of M07 include the neu- 
trino micro-physics necessary for direct application to PNSs. 
Their calculations were, however, restricted to monopolar 
field structu re. A complementary set of calculations was car- 
ried out by iBucciantini et all ( 20061 ) (hereafter B06) , who 
studied aligned dipolar (axisymmetric) non-relativistic and 
relativistic MHD winds from neutron stars assuming an adi- 
abatic equation of state. M07 showed that their results could 
be combined with those of B06 to provide evolutionary mod- 
els for early magnetar spin-down, including the transition 
from non-relativistic to relativistic outflows as the neutrino- 
driven mass loss decreases. 

Figure [1] shows the energy loss rate E and magnetiza- 
tion olc as a function of time from the evolutionary cal- 
culations of M07 for several combinations of (dipole) mag- 
netic field strengths and magnetar birth period. The values 



of B w 10 1 



10 G and P w 1 — 2 ms are chosen to be 



characteristic of PNSs that produce conditions suitable for 
producing GRBs or hyper-energetic supernovae. The mag- 
netization in Figure [T] is defined by 

_3 



(1) 



where $b is the total open magnetic flux per 4-7T steradian 
(Michel 1969), M is the mass loss rate, and olc is evalu- 
ated at the light cylinder. Winds with olc < 1 are non- 
relativistic while those with olc > 1 are relativistic. The 
calculations in Figure [1] assume that the PNS luminosity 
decreases in tim e oc t~ unti l t = 4 s, motivated by the cool- 
ing evolution of IPons et alJ (jl999l ). We note, however, that 
the cooling of strongly magnetized rapidly rotating PNSs is 
not well understood, which is a source of uncertainty in the 
determination of M(t), E(t) and OLc(t)- 

The energy loss rates E in Figure [T] are a factor of ~ 10 
larger than the "vacuum dipole" rate for the given value of f2 
and B. There are two reasons for this. First, when the wind 
is non-relativistic (olc <J 1), the energy loss rate is larger 
than in the relativistic limit, with E oc M 1/3 (M 3/5 ) for non- 
relativistic magnetized win ds with a monopole (di pole) mag- 
netic field structure (e.g.. [Thompson et al.l [20041 ). In addi- 
tion, the large mass loss rates accompanying PNS spin-down 
lead to excess open magnetic flux which enhances the spin- 
down of the PNS (see the simulations of B06). This is true 
even when the wind is moderately relativistic (olc > 1). 
The large energy loss rates shown in Figure[T]are sufficient to 
lead to appreciable spin-down of the PNS during the Kelvin- 
Helmholtz epoch. For example, for the model with P — 1 ms 
and B = 3 x 10 15 G in Figure [l] (solid line), the PNS loses 
~ 80% of its rotational energy in the first 40 seconds. This 
efficient spin-down is largely responsible for the fact that E 
decreases in time as the PNS cools (see Figure [1} 



2 Two additional effects contribute to the decline in E with 
time. First, as the PNS cools, the mass loss rate M decreases. 
In the non-relativistic limit, the energy loss rate is proportional 
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As the PNS cools, the neutrino-driven mass loss de- 
creases in time. This in turn causes a transition from a non- 
relativistic to relativistic wind, as shown explicitly in the 
plot of ULc(t) in Figure [1] These calculations of glc are 
based on equatorial spin-down models (M07), which prob- 
ably underestimate the angle- averaged a in the wind by a 
factor of few (B06). Nonetheless, the evolution from a mod- 
erately mass-loaded marginally relativistic wind (glc ~ 1) 
to a highly magnetized Poynting flux dominated outflow 
(glc S> 1) is expected to be generic for cooling magnetars. 

As we show in the next section, the impact of the mag- 
netar on its host star depends critically on the strength of 
the magnetic field in the bubble created by the magnetar 
wind; the generation rate of the field in the bubble is in 
turn determined by the magnetization a of the wind at large 
radii. In non-relativistic winds, the magnetic energy and ki- 
netic energy are in approx imate equipartition at larg e radii, 



with Eraag « 2Eki n (e.g.. lLamers fc Cassinellil ^999). One- 
dimensional models of ideal relativistic winds, however, find 

1/3 

that the asymptotic Lorentz factor of the wind is as g lc 
and the asymptotic magn etization is a « ct^q (|Michelll 19691 . 
iGoldreich fc Julian! 11970 ) so that most of the energy re- 
mains in the magnetic field at large radii. These results ap- 
ply in the limit of glc 1. Relativistic MHD simulations 
l|Bucciantini et all 120071 ') show that for intermediate values 
of glc < 20, a reasonable fraction of the magnetic energy is 
converted into kinetic energy at large distances, with rough 
equipartition obtaining by ~ 10 4 stellar radii. 

In the limit of very high glc, studies of PWNe (e.g. the 
Crab Pulsar) find that the wind must have low g ~ 10~ 2 
at lar ge radii fe.g.. lKennel fc Coroniti|[l984l . iBegelman fc Lil 
1992). Although there is no consensus on the mechanism 
responsible for the inferred decrease in pulsar wind mag- 
netization at large radii, a prominent class of models relies 
on magnetic dissipation in the relativistic outflow over a 
large radial distance (e.g., Cor oniti 199Cj|; lLvubarskv fc Kirkl 



l200ll : iKirk fc Skjaeraasenll2003l ). The physical conditions in 
proto-magnetar winds are quite different from those in pul- 
sar winds (e.g., they are much denser so that there is no 
charge starvation). In addition, the distance to the termina- 
tion shock is much smaller in the SN confined winds from 
young magnetars, ~ 10 light cylinder radii (see below) com- 
pared to more than 10 4 light cylinder radii in PWNe and 
in pulsar-Be star binaries. The reduced flow time between 
the light cylinder and the termination shock may mean that 
dissipation of magnetic energy in young magnetar winds is 
less complete than in pulsar winds. As a result, we suspect 
that the rate of injection of magnetic energy into bubbles 
created by protomagnetars may be significantly larger than 
that inferred in the PWNe context. Given the uncertain- 
ties, however, we treat the magnetization in the outflow, 
expressed as the ratio of the magnetic energy injection to 
the total power (E mag / Etot), as a free parameter in this pa- 
per, bearing in mind the generic evolution from glc ~ 1 to 
glc > 1 in Figure Q] 



to M 3 0,6 and thus decreases as well (this is relevant for the 
P = 1 ms, B = 3 X 10 15 (solid line) and P = 1 ms, B = 10 15 G 
(dot-dashed) models in Figure [l] at early times). The decreasing 
mass loss rate also decreases the fraction of open magnetic flux 
and thus E. 



The models shown in Figure \T\ assume that the wind 
from the central magnetar is freely expanding into a cav- 
ity evacuated by the outgoing SN shock. Formally, this re- 
quires that the radius of the fast magnetosonic point must 
be smaller than the radius of the SN shock; the latter is 
R s ~ 10 9 cm in the first few seconds, which is indeed larger 
than the typical distance to the fast surface of ~ 10 — 40 
neutron star radii (B06, for a millisecond rotator). As the 
freely expanding wind moves out, it interacts with the sur- 
rounding SN shock and previously shocked wind material. 
More precisely, the wind will reach a termination shock at 
which its kinetic energy is thermalized and the magnetic 
field is compressed. A correct determination of the size of the 
termination shock requi res a full MHD model of the wind- 
bubble interaction ('e.g.. lDel Zanna et al.|[2004T ). As a rough 
guide to the relevant scales, however, we note that in the 
simple case of a constant M and E, glc ~ 1 wind moving 
into a spherically symmetric bubble, the termination shock 
is located at a radius R t ~ RsiRs/ct) 1 ^ 2 ~ 0.1-R 3 ~ 10 s cm 
where t is the time since the onset of the wind (in sec). For 
Rt < R < R s , the wind develops into a bubble of plasma 
and magnetic field confined by the SN shock and host star. 

2.2 The Bubble Structure 

If one neglects plasma flow inside the bubble, a simple solu- 
tion for the structure inside the bubble (R t < R < R s ) can 
be obtained in the case of the predominantly toroidal mag- 
netic field ex pected at large radii in the wind. This solution 
was found bv lBegelman fc Lil (|l992T) . We reproduce several 
of its fea tures here because they are important to our model. 

The lBegelman fc Lil l|l992T ) solution will be valid as long 
as typical flow speeds do not exceed the local sound speed. 
In the case of a relativistic magnetized bubble the sound 
speed ranges from c/y3 to c. It is possible that, close to 
the termin ation shock, post shock flow can move with high 
velocities (|Del Zanna et all l2004h , but in the bulk of the 
bubble, typical speeds are expected to be a small fraction of 
c, unless the cavity itself expands at a considerable fraction 
of the speed of light. Indeed, as long as the expansion veloc- 
ity of the shell is small compared to the sound speed inside 
the bubble, the plasma inside will always relax to pressure 
equilibrium, independent of the energy distribution in the 
wind (be it primarily polar as for a non-relativistic wind 
or primarily equatorial as for a relativistic wind). Neglect- 
ing the fluid flow, the structure is given by the assumption 
of magnetohydrostatic equilibrium. Assuming axisymmetry, 
the momentum equations become: 



d_ 

dz 



p + 



8tt 



0. 



d_ 

dr 



P + 



8tt 



+ 



47rr 



0. 



(2) 



where r is the cylindrical radius, p is the pressure, and B 
the toroidal magnetic field in the bubble. The first equation 
simply states that isobaric surfaces are coaxial cylinders. 

If entropy is constant along each flow streamline in the 
bubble then the continuity equation can be written as: 



® 1 i/r \ . d , i/r % „ 
— (p rvr) + ^-(p rv z ) = 0. 
or oz 



(3) 



where F is the adiabatic index of the fluid. Comparing this 
with the flux-freezing condition for the toroidal magnetic 
field yields 
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(4) 



For the case of a relativistic plasma (F = 4/3), equation Q 
can be used in the r-momentum equation to find 



Pn 



B 2 



9p n r 2 



C 2 ' 8tt 16C 3 # 2 ' 
where £ is the solution of the following equation: 



« + 



32H 



:) 2 -C 3 



0. 



(5) 



(6) 



The solution for the pressure in the bubble given by 
equations ([2j)- ((6j) depends on two parameters. One of these, 
the pressure on the axis p„, determines the overall mag- 
nitude of the pressure in the bubble. The other, the scale 
height H of the pressure distribution, determines the pres- 
sure stratification in the bubble. In Figure [2] we plot the 
normalized pressure profile derived from the solution of the 
above equations. The total pressure is higher along the axis 
(r = 0) and asymptotically decreases as r~ 2 . The region 
close to the axis contains a low a plasma and is essentially 
pressure dominated, while at larger distances the plasma is 
magnetically dominated, and the ratio of magnetic to ther- 
mal pressure increases linearly with the distance. Equipar- 
tition is reached for r/H ~ 2. The results in Figure [2] as- 
sume a relativistic plasma with F = 4/3, which corresponds 
to olc > 1 in Figure [T] The magnetar wind may be non- 
relativistic at very early times, so that V — 5/3 is more 
appropriate. For T = 5/3 the pressure profiles are qualita- 
tively similar to those in Figure [5] although the scale height 
H is a factor of w 2 smaller for a given ratio of magnetic 
to total energy in the bubble. For simplicity, we simply set 
r = 4/3 in all of our calculations. 

The scale height H and the asymmetry of the pressure 
distribution can be expressed in terms of the ratio of the 
magnetic energy to total energy in the bubble. To quantify 
this effect, consider a spherical bubble of radius R and total 
energy E. The pressure along the axis is given by 



E 



10 51 ergs 



)( 



R 



10 9 cm 



ergs cm 



(7) 



The dimensionless number P is the pressure on the axis 
relative to that in an unmagnetized bubble. Figure [3] shows 
P and the scale height H/R as a function of E mag / Etot, the 
ratio of the magnetic to total energy in the bubble (similar 
results are obtained for the self-similar solutions described 
in the Appendix; see, e.g., Figure [A"2)l . Magnetized bubbles 
have P 3> 1 and H <C R (where R is the radius of the 
bubble, not the cylindrical radius within the bubble used 
above and in Figure HJ. Figure [2] shows that, due to the 
pinching effect of the toroidal magnetic field, the pressure in 
the bubble in this case will be concentrated along the axis 
and so the bubble will expand asymmetrically. By contrast, 
very weakly magnetized bubbles have H > R and roughly 
constant pressure throughout. Note that a magnetization 
of E mag I Etot J> 0.1 is required to make H <C R and the 
pressure distribution in the bubble relatively asymmetric. 

We now calculate how the swept-up shell in the host 
star responds to the pressure produced by the magnetized 
bubble created by the central magnetar. 
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Figure 2. Pressure structure for a relativistic plasma containing 
to roidal magnetic field in magneto- hydrostatic equilibrium (based 
on lBegelman fc Lilll992h . The solid line (A) is the total pressure, 
normalized to the value on the axis. The dotted line (B) is the 
magnetic pressure, normalized to the total pressure on the axis. 
The dashed line (C) is the ratio of the magnetic to the total pres- 
sure, while the dash-dotted line (D) is the ratio of the magnetic 
to the thermal pressure. Bubbles with weak magnetic fields have 
large values of H relative to the size of the bubble (see Figure 
(3)l and thus only the r H part of this plot is applicable: as a 
result the pressure is relative uniform and the system will expand 
spherically. By contrast, bubbles with appreciable magnetic fields 
have smaller values of H and thus the pressure on the axis is sig- 
nificantly larger than the pressure near the equator. Such bubbles 
will expand asymmetrically. 




0.6 



Figure 3. Dimensionless pressure P (see eq. [7]) and scale-height 
H as a function of the magnetic energy in the bubble, for the 
case of spherical bubble. P is the ratio of the pressure on axis 
to the pressure in an unmagnetized spherical bubble. For even 
modest magnetic energy, the pressure distribution becomes highly 
anisotropic with P>1 and H R. 
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3 THE THIN-SHELL EQUATIONS 

The equations describing the evolution of an axisymmet- 
ric bubble can be greatly simplified if one assumes that the 
swept-up material is confined in a thin-shell, so that the dy- 
namics can be described in the "thin-shell approximation." 
This approximation is valid as long as the thickness of the 
shell is small compared to the radius of curvature. The thin- 
shell equations account for conservation of mass and momen- 
tum. A detailed derivation of the equations can be found in 
iGiulianil l| 19821 ) where corrections for the thickness of the 
shell are also discussed. In the case of infinitesimally thin 
shell they are given by: 

tan£ = -- — , (8) 



u± 
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where £ is the angle between the radial direction and the 
normal to the shell surface, A is the effective area of each 
element of the shell, and a is the surface density. The suffixes 
_L and |] represent the directions perpendicular and parallel 
to the shell (and they are time dependent), while i and o 
label the conditions at the inner and outer edge of the shell. 
The velocity ui = u± is the expansion velocity of the shell 
perpendicular to itself, u« is the expansion velocity parallel 
to itself, and v\\ is the flow velocity along the shell. 

Equation (|ll|l represents conservation of mass along the 
shell, while equation (|12|l describes momentum conservation 
in the shell. Both equations include a flux of the relevant 
quantity along the shell itself, and so urce terms due to the 
inner and outer media. As discussed in IGiulianil l| 1982b . these 
equations employ a Lagrangian remapping along the shell, 
and can be applied only as long as R(9) is a single valued 
function of 6. 

The evolution of the thin shell depends on the force 
driving it ( "i" ) and on the conditions in the external medium 
("o"). The solution in H2.2I describes the inner conditions 
used in this study. In many cases of interest, the outer 
medium is sufficiently cold and the magnetic field is suffi- 
ciently weak that their contribution to the pressure term 
in equation (|12|l can be neglected. In our case, the outer 
medium is the outer part of the progenitor star from which 
the magnetar formed (see Q; we do not neglect the thermal 
pressure, but we do assume that the progenitor is unmagne- 
tized. 

Given the evolution of H and p„ with time (calculated 
below), equations (f8)l- (|13[l were solved under the assumption 
of axisymmetry, to determine the evolution of the shape of 



the shell with time. We were not able to cast these equa- 
tions in full upwind form, because of the presence of a term 
describing the advection of the shell curvature, which is not 
constant and changes in time. This requires adding some ar- 
tificial viscosity in order to avoid the numerical growth of 
perturbations. 

One cautionary comment about the thin-shell model is 
in order. Equations (f8)l- (|13[l are momentum conserving, not 
energy conserving, in the sense that a shell expanding into 
an ambient medium has constant momentum and thus its 
energy decreases in time. The equations do conserve energy, 
however, in the sense that the work done by the interior 
bubble is self-consistently supplied to the shell (see §[3}, but 
some of this energy is then lost as the shell expands and 
sweeps out into the ambient medium. Unfortunately, it is 
not possible to conserve both momentum and energy in the 
time dependent, thin-shell approximation (by contrast, in 
the self similar case discussed in the Appendix, one can sat- 
isfy both requirements, but in this case the time evolution 
is factored out of the equations). One consequence of this is 
that the calculations that follow probably evolve somewhat 
more slowly than would a true SN shock expanding into its 
host star, although we are confident that our conclusions 
about generating asymmetric bubbles are robust (e.g., the 
self-similar solutions in the Appendix show similar asymme- 
try). 



4 A MAGNETAR IN A BOTTLE 

In this section we use the thin-shell model to calculate the 
evolution of the magnetized bubble inflated by a central 
magnetar. As was pointed out in §2.1, one needs to derive 
the internal pressure distribution in the bubble in order to 
solve for the dynamics of the shell. In particular, one needs 
to know the value p n of the total pressure on the axis and 
the value of the scale height H of the pressure distribution. 
Once these two parameters are known it is possible to derive 
the pressure throughout the bubble, in particular its value 
at the inner edge of the shell. One can show that given the 
shape of the shell bounding the bubble, p n , and H , the total 
energy E to t, the magnetic energy E mag , and the magnetic 
flux <E> inside the bubble itself are uniquely defined, where 



Etot — / p„J r (r/H,z)dv, 
Jv 

E m ag = / P,iG (r / H , z)dv , 

Jv 



$ = 



/p^Q(r/H,z)da, 



(14) 
(15) 
(16) 



and where V is the volume of the bubble and A is the area 
in the r — z plane, delimited by the shell. The dimensionless 
functions T,Q,Q are given in terms of cylindrical coordi- 
nates, and can be derived from the pressure and magnetic 
field given by equations (E}-© in § 2.1. 

In order to compute the evolution of the internal struc- 
ture in the bubble we subdivided each time step (dt) of the 
shell evolution into two sub-steps. In the first sub-step, that 
we call adiabatic, we neglect injection of energy and mag- 
netic field by the central source, and we compute the adia- 
batic losses due to expansion according to: 
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dE tot = J P dV, (17) 

where p is the total pressure along the shell surface S and 
dV is the volume increment that results from the evolution 
of the shell surface. Once the adiabatic losses are known 
one can derive the new value for the total energy in the 
bubble. During this adiabatic step the magnetic flux remains 
constant. After the adiabatic step, the new values of p n and 
H are re-evaluated by solving the following equations: 

E to t,a = E tot - dE to t = J p n F(r/H, z)dv, (18) 
* = / ^/p7 l Q{r/H,z)da, (19) 

J A 

where the integrals are computed using the values of V and 
A after the expansion. Using the new values of p n and H, we 
need to recompute the new magnetic energy inside the bub- 
ble E„ la g,a, because adiabatic losses act on the total energy. 
This is done using equation (|15[) . 

In the second sub-step, that we call the injection step, 
the shape of the bubble is assumed to be fixed and we com- 
pute the new values of the total energy and the magnetic 
energy given the rate of total energy and magnetic energy 
injection by the central magnetar. The two equations to be 
solved for p n and H are: 

E t ot,a + Etotdt = / p n T{r/H,z)dv, (20) 
Jv 

E mag ,a + E mag dt = / p n Q (r I H ', z)dv , (21) 
Jv 

and once p n and H are known we can also recompute the 
magnetic flux $, which will be needed in the next time step. 
With this method we determine the evolution of the pressure 
on the inner edge of the shell as a function of time given 
E t ot{t) and E mag (t) (by, e.g., the results of Figure [T|. 

Based on modeling the spectra of supernovae associated 
with nearby GRBs, there are some indications that GRBs 
arise from v ery massive stars wit h ZAMS masses of M ~ 
40Mq (e.g., Ilwamoto et al.|[l998f ). There are also observa- 
tional indications that Ga lactic magnetars are formed from 
comp arably massive stars i|Gaensler et all 2005 , Muno et al.l 
2006). We thus consider the evolution of a magnetized bub- 
ble inside a p rogenitor star o f 35 Mp, using the progeni- 
tor models of IWooslev et al.1 i|2002l ). We have also consid- 
ered lower progenitor masses down to « 11 A/q, which may 
be more appropriat e for the progenitors of X-ray flashes 
l|Mazzali et alJuOOfJ ). We find little difference in the results 
for different progenitors, at the level of quantitative detail to 
which we aspire. The most significant effect is that for fixed 
energy injection, the bubble expands more slowly for more 
massive progenitors. The relative insensitivity to progenitor 
mass can in part can be understood by noting that the self- 
similar solutions described in the Appendix show explicitly 
that the elongation of the bubble depends only weakly on 
the density profile of the ambient medium. 

As discussed in §2, our model of the magnetar wind as- 
sumes that it is expanding into a cavity evacuated by the 
outgoing SN shock. To initialize our simulations, we thus 
carve out a spherical cavity with a radius of 10 9 cm inside 
our progenitor, corresponding to the region of infall in the 
first ~ 1 sec. We assume that this cavity is bounded by a 



thin shell whose mass is equal to the mass that originally 
was in the cavity region minus IAMq (the canonical mass 
for a neutron star). In all of our simulations, time is defined 
after core bounce and the simulation starts 1 second after 
core bounce. Moreover we impart to the shell an outward 
velocity so that the total shell energy at the beginning is 
10 51 ergs, enough to trigger a SN. If instead one assumes 
an initially stationary shell, the evolution is essentially un- 
changed for weakly magnetized bubbles because the pressure 
of the bubble is relatively isotropic (this assumes that the 
magnetar wind extracts at least ~ 10 51 ergs at early times, 
as is the case in the models shown in Figure [TJ. For strong 
magnetization, the elongation of the bubble along the axis 
is also nearly independent of the initial shell energy. How- 
ever, for large E mag /Etot, the pressure in the bubble near 
the equator can be so small that infall cannot be prevented. 
To model this case, a full hydrodynamic solution is required. 

We follow the evolution of the shell and interior bub- 
ble to large distances, into the hydrogen envelope of the 
progenitor. For GRB progenitors, the hot plasma confined 
inside will emerge into the circumstellar medium once the 
shell surface reaches the outer edge of the helium core. The 
initial material that emerges will probably only move with 
modest Lorentz factor. Subsequent material will, however, 
rapidly accelerate through the channel carved by the mag- 
netar wind, reaching asymptotic Lorentz factors set roughly 
by the enthalpy of the material in the bubble (assuming 
that E mag <J Ethermai in the bubble). This phase of evolu- 
tion cannot be studied using the thin shell-approximation, 
but requires full relativistic MHD simulations. Nonetheless, 
it appears natural that a highly relativistic and collimated 
outflow will emerge out of the cavity carved by the early 
magnetized bubble. 

4.1 Results 

In Figure [3] we show the results of a series of simulations 
for different fixed values of Emag/Etot, the ratio of the 
Poynting flux to the total power injected by the magne- 
tar. In all of the calculations in Figure 3J we assume that 
the total power supplied by the central source is given by 
Etot = 10 51 (t/1 s) _1 ergs -1 , which is a reasonable approxi- 
mation to the lower power solutions in Figure [T] (note that 
we neglect the possibility of early injection of energy and 
start our simulation 1 second after core bounce). Note that 
in this case equal energy is supplied per decade in time. 

For values of Emag/Etot ~ 0.1, the pressure distribution 
inside the bubble is relatively spherical (Figs. and so 

the surrounding shell becomes only modestly asymmetric. 
Most of the energy supplied by the central magnetar in this 
case is transferred to the surrounding SN shock and host 
star. Low magnetization bubbles of this kind would thus 
likely produce a mildly asymmetric hyper-energetic SNe, but 
it appears unlikely that the relativistic material supplied by 
the magnetar can easily escape its host star. For larger values 
of Emag/Etot, the shell evolves more asymmetrically because 
most of the pressure is exerted along the axis for magnetized 
bubbles. By Emag/Etot > 0.3, there is clear evidence for a 
very elongated channel driven through the host star by the 
anisotropic pressure of the central bubble. The shell reaches 
the outer edge of the progenitor (~ 2 ■ 10 10 cm) after ss 
5—10 sec. At this point the ambient density drops to typical 
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Figure 4. Evolution of a magnetized bubble inside a 35Mq progenitor for a central source with Etot = 10 51 (t/1 s) — 1 ergss -1 . The 
initial time is Is; the contours describe the shape of the shell at Is time intervals for the first three panels and 0.5s time intervals for the 
last. From left to right, the solutions correspond to increasing the magnetization of the bubble, with E mag / Etot = 0.1,0.2,0.3, &0.4, 
respectively. For Emag / Etot > 0.3 most of the pressure of the inner bubble is exerted close to the axis (see Figs. [2] & [3j , which leads to 
the asymmetric evolution of the bounding shell. 



circumstellar values, and the shell will rapidly blow out of 
the star. The highly relativistic material contained in the 
interior bubble can now flow relatively unimpeded out of the 
host star, forming a relativistic jet; it is natural to associate 
these models with the production of a GRB. 

Figure [5] shows the evolution of the thin shell for a 
more energetic, but more rapidly decaying, central source 
with Etot = 10 52 {t/1 s) -2 ergs -1 , which is an approxima- 
tion to the higher power solutions in Figure [1] We consider 
E mag /E t ot = 0.2 (left) and E mag /E to t = 0.3 (right). Note 
that in this case, most of the energy is supplied to the bub- 
ble at early times and so the evolution of the system is sim- 
ilar to the case of a magnetic bomb with a fixed energy 
of ~ 10 52 ergs in the bubble. The evolution of the shell in 
Figure [5] is qualitatively similar to that of the lower power 
solutions shown in Figure 3J although the bubble evolves 
more rapidly because of the more energetic central source. 
One consequence of this more rapid evolution is that the 
shell velocity is closer to c, implying that the assumption of 
magneto-hydrostatic equilibrium used to derive the interior 
structure is likely to be less accurate than in the case of the 
weaker power sources in Figure [3] 

For PNSs with rotation periods longer than the val- 
ues of ~ 1 — 2 ms considered in Figure [T] the energy injec- 
tion rate will be lower and approximately constant at early 
times because the spindown time is longer than the Kelvin- 
Helmholz time of w 10 — 100 s. To investigate this limit, we 
considered the evolution of a bubble with a constant energy 
injection rate of Etot ~ 10 50 erg/s. Elongation analogous 
to that shown in Figures [3] & [S] can be achieved, although 
somewhat higher magnetization is required. An asymmetry 
similar to the E mag / E to t = 0.2 solution in Figure [4] requires 
E ma g/Etot = 0.3 for this lower Etot and takes a somewhat 
longer time ~ 20 sec to develop. This example highlights 
that lower power sources - which can originate from more 
modestly rotating PNSs - can still lead to asymmetric bub- 
bles because the energy per unit solid angle along the pole is 
significant even for modest E to t ~ 10 49 - 10 50 ergss . Such 
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Figure 5. Evolution of a magnetized bubble inside a 35Af0 
progenitor for a central source with a total spindown power of 
10 52 (i/ls) -2 erg s — 1 ; the initial time is Is. Contours represent 
the shape of the shell at Is intervals. Left: Emag/Etot = 0.2 Right: 
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: 0.3. 
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sources may form asymmetric SN and, in some cases, 
long-duration GRBs or X-ray flashes. 

An approximate analytic understanding of the late-time 
structure of the shell shown in Figures [4] & [5] can be ob- 
tained by assuming that most of the energy is released in- 
side the bubble before it has time to expand significantly - so 
that its shape is still approximately spherical and adiabatic 
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losses can be neglected. In this case most of the pressure 
will reside in a region along the axis whose opening angle 
is ~ H/R and the shell will expand rapidly along the axis 
within this solid angle. Figure [3] gives the relationship be- 
tween H/R and the magnetization of the bubble, which can 
be used to estimate the opening angle of the resulting "jet" 
at late-times. For example, Figure [3] shows that H/R ~ 0.1 
for E m ag / Etot ~ 0.3 — 0.4, which is reasonably consistent 
with the angular collimation of ~ 0.1 rad in Figure 2] It is 
also worth noting that the high axial pressure produced by 
a magnetized bubble with energy E leads to an expansion 
along the axis that is quantitatively similar to that produced 
by a jet with kinetic power 

where we have assumed that the angular collimation 8 ~ 
H/R and that the dimensionless axial pressure P « R/H 
(which is reasonable for H > 0.03 J?; Figure [3J. 

The results in Figures U & [S] assume that Emag/Etot 
is independent of time. This may be a poor approximation 
given the strong evolution in the magnetization of the wind 
as a function of time at small radii (Fig. [TJ. Indeed, one 
might naively expect that Emag/Etot would increase in time 
on a timescale of a few sec, given the evolution of tJLc(t) 
for magnetar winds. If this is correct, the magnetar-created 
bubble may initially impart its energy relatively spherically, 
enhancing the energy of the outgoing SN shock (as in the 
left panel of Figure [4|. At late times, however, the bub- 
ble will drive a jet out along the axis (as in the right-most 
panel in Figure 3]). To explore this scenario, we carried out 
a series of simulations starting with a small initial value of 
a ss 0.05—0.1 and then increasing a in time as <jlc increases, 
up to the equipartition value of a « 0.5. As expected, in the 
first few seconds the evolution of the shell was quite spher- 
ical, closely resembling the a = 0.1 panel of Fig. [4] while 
at late times the shell evolves into a more elongated struc- 
ture analogous to the a = 0.3 — 0.4 panels of Fig. [4] In this 
scenario, the different panels in Fig. [4] qualitatively describe 
different snapshots in time for the evolution of a shell driven 
by a wind with increasing magnetization. This suggests that 
the increasing magnetization of the magnetar wind provides 
a way of tapping the spindown energy to generate a mildly 
asymmetric hyper-energetic SN, while at the same time cre- 
ating an axial cavity along which relativistic material can 
escape, powering a GRB. 

Throughout this paper, we have considered only models 
in which the progenitor has a spherical density profile, in or- 
der to understand the importance of magnetic stresses on the 
asymmetric evolution of the stellar wind bubble. However 
for rotating stars the density in the polar region is expected 
to be smaller than at the equator. This can also facilitate 
collimation. A full investigation of the combined collimation 
induced by magnetic stresses and a non-spherically symmet- 
ric stellar profile is beyond the scope of this paper. We have, 
however, carried out a few preliminary calculations investi- 
gating the two extremes: a low density polar channel whose 
opening angle is either smaller or bigger than H/R vo i ar . Our 
results show that at low magnetization the presence of a low 
density channel along the rotation axis can lead to significant 
collimation (well in excess of the mild asphericity shown in 



Figure [4] at low magnetization) , while for higher magnetiza- 
tion the collimation is primarily due to the magnetic stresses 
we have focused on in this paper. 



5 DISCUSSION 

In this paper we have calculated the evolution of a magne- 
tized bubble formed inside an exploding massive star. Our 
motivation is to understand the impact of a newly born mil- 
lisecond magnetar on its surrounding stellar envelope, and in 
particular to determine the conditions under which magne- 
tar birth can produce the collimated outflows inferred from 
observations of long-duration GRBs. 

Neutron stars born with P ~ 1 ms and B ~ 10 15 - 10 16 
G lose ~ 10 51 - 10 52 ergs in ~ 10 — 100 sec in a mag- 
netized wind that becomes increasingly relativistic at late 
times (Figure [1]). This energy forms a bubble of plasma and 
magnetic field confined by the inertia of the surrounding 
star. If the material supplied to the bubble has low mag- 
netization, Em a g/Etot 0.1, the bubble expands relatively 
spherically (Figure [3} and most of the energy of the spin- 
ning down neutron star is transferred to the surrounding SN 
shock, plausibly powering a hyper-energetic SN, but proba- 
bly not a GRB. By contrast, for more appreciable magne- 
tization, Emag/Etot }> 0.3, the bubble created by the mag- 
netar rapidly becomes highly elongated along the rotation 
axis of the system as a result of the internal pressure dis- 
tribution (§2.1), forming a cavity out of which the late-time 
relativistic wind from the magnetar can escape (Figure [4] & 
0). We suggest that this is plausible physical mechanism for 
forming collimated relativistic outflows from newly formed 
millisecond magnetars(f] 

This mechanism works even if the host star is spheri- 
cally symmetric. In addition, even if most of the wind energy 
flux is concentrated in the equatorial plane (as is expected 
for relativistic winds from neutron stars), the surrounding 
bubble will still reach magnetohydrostatic equilibrium and 
will thus elongate along the axis due to magnetic stresses as 
we have described. Finally, we note that it is not necessary 
to wait until late times, when the magnetar wind is rela- 
tivistic, in order for the bubble to evolve highly asymmetri- 
cally. Collimation can be achieved in the early mass loaded 
phase, provided only that the bubble contains a sufficient 
toroidal magnetic field. This mechanism may then operate 
in additio n to hoop-stress collimation of the non-relativistic 
flow (e.g., iBucciantini et all 120061 : lUzdenskv fc MacFadvenl 
2006). This early time non-relativistic phase cannot by it- 
self produce a GRB, but can create a channel out of which 
the later relativistic wind emerges. Such a channel might 
also provide boundary conditions conducive to the acceler- 
ation of the wind and the conversion of magnetic energy 

3 Although we have framed much of this paper in terms of the 
magnetar model for GRBs, magnetized outflows from an accre- 
tion disk around a neutron star or black hole would produce asym- 
metric magnetized bubbles akin to those considered in this paper. 
Whether they would play a similar role in driving collimated flows 
inside a supernova depends on the details of the disks' mass and 
magnetic energy losses as a function of time, as well as uncer- 
tainties in the disk winds' a similar to those encountered in the 
neutron star winds considered here. 
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into kinetic energy (|Komissarov et all I2007T ) . Our calcula- 
tions show that for the expected magnetar energy loss rates, 
a collimated cavity is formed after ~ 10 sec (Fig. [4} . At this 
point, magnetar winds have olc ~ 100 (Fig.[l|, in the range 
required to account for GRBs. 

Becaus e the birth rate of magne tars (~ 10% of neutron 
stars; e.g., iKouveliotou et ail Il994l ) is significantly larger 
tha n the rate of GRBs (~ 0.1 — 1% of massive stellar deaths; 
e.g.. [Podsiadlo wski et al.1l2004 ). most magnetar births can- 
not produce standard long-duration GRBs. This is presum- 
ably either because an extended stellar envelope inhibits the 
escape of a collimated outflow or because most magnetars 
are born rotating more slowly than the millisecond rota- 
tors we have focused on in this paper. For more modestly 
rotating PNSs, the asymmetric expansion of a magnetized 
bubble could contribute to the inferred asymmetr y of many 
core-collapse supernovae (e.g., IWang et al.ll200lT ). In addi- 
tion, a PNS with, e.g., P « 4 ms and B » 3 x 10 15 G 
has a rotational energy of « 10 ergs and a spindown time 
of ~ 1 day. The birth of such a neutron star would not 
produce a hyper-energetic SN or a canonical GRB. How- 
ever, if the bubble created by the magnetar is sufficiently 
magnetized, it would evolve asymmetrically in a manner 
similar to the calculations shown in Figures |4] & This 
could produce a long-duration transient analo gous to the X- 
ray f l ash 060218 associated with SN 2006aj |Mazzali et all 
2006); ISoderberg et all [2006; we should note, however, that 
many X-ray flashes may have lower inferred energies be- 
cause of viewing angle effe cts rather than bei ng intrinsically 
less energetic events; e.g., (|Granot et aUl2"005h . The remnant 
Cass A, with its strong jet/counter-jet morphology (e.g., 
I Hwang et al. I [20041 ). may be an example of an asymmetric 
explosion driven and shaped by a magnetized wind a ccom - 
panying magnetar birth. Indeed, IChakrabartv et all (1200 ll ) 
suggested that the central X-ray point source in Cass A is a 
magnetar. 

The thin-shell calculations described in this paper as- 
sume that the magnetar wind expands into an initially 
spherical cavity created by an outgoing SN shock. This re- 
quires that the spindown time of the magnetar is at least 
somewhat longer than the time required to initiate the stel- 
lar explosion (i.e., > 1 — 2 sec). Our assumption of a "suc- 
cessful" SN explosion does not, of course, preclude that the 
explosion itself is magneto-centrifugally driven, as in the 
force-free model for the collimated explosion of a star by a 
newly-formed magneta r in an otherwise "failed" SN (e.g ., 
lOstriker fc Gunnl Il97ll or lUzdenskv fc MacFadvenl 120071 ). 
However, one interesting problem not addressed by our cal- 
culations is the spindown of the magnetar and the evolu- 
tion of its surrounding bubble if the initial explosion is pri- 
maril y bi polar (see, e.g., the s imulations of Moi seenko et al.l 
120061 and lBurrows et al"1 [2007). Late-time collimation of rel- 
ativistic material in this context may be modified by the 
large inertia of the accreting stellar envelope (or fallback 
material) in the equator of the star (see al so the related ar- 
guments of lUzdenskv fc MacFadvenl 120071 ). In addition, it 
is worth noting that if the outflow always has high mag- 
netization, our calculations suggest that, because most of 
the pressure will be exerted along the axis, there could be 
a collimated GRB but no associated equatorial explosion. 
This could account for the re cently discovered supernova- 
less GRBs (|Fvnbo et alj|2006h . 



One of the uncertainties associated with our calcula- 
tions is that the magnetization of the material supplied to 
the surrounding bubble is difficult to calculate. Magnetic 
energy has to be supplied to the bubble relatively rapidly, 
with a = Emag/Etot J> 0.2 at the termination shock in our 
models that show significant collimation. Observations of 
PWNe suggest quite low a ~ 0.01 at the termination shock, 
which would imply that there is insufficient time to build up 
the anisotropic magnetic stress needed to drive aspherical 
expansion of the surrounding stellar envelope. However, we 
suspect that the confined bubbles around newly formed mag- 
netars will have higher magnetization at their termination 
shocks than has been inferred in PWNe and in pulsar-Be 
star binaries. This is because the distance to the termina- 
tion shock is only ~ 10 light cylinder radii in our problem, 
relative to > 10 4 light cylinder radii in the systems where 
we have direct observational constraints. As a result, there 
is less time for the magnetic field in the wind to dissipate, 
plausibly leading to higher magnetization. 

All of the calculations described in this paper are based 
on the thin-shell approximation. This model is useful for 
demonstrating that magnetar birth can produce conditions 
conducive to the formation of a collimated outflow that can 
emerge out of the host star. However, modeling this pro- 
cess in detail is beyond the scope of the present simplified 
calculations, and will require full relativistic MHD simula- 
tions. Indeed, it is our intention to use the results of the 
present paper as a guide for more realistic simulations. Such 
calculations are necessary to determine the fraction of the 
spindown energy that goes into a relatively spherical explo- 
sion of the host star relative to the energy that flows out 
of the collimated cavity. Quantifying this is important for 
understanding the conditions under which magnetar birth 
might produce b oth a hyper-energetic SN and a GRB, as 
is observed fe.g.. IWooslev fc Blooml 120061 ) . We have specu- 
lated in §4 that this can occur if the surrounding bubble be- 
comes progressively more magnetized as the magnetar spins 
down, but multi-dimensional simulations are needed to as- 
sess this (and to understand if it is a necessary condition). 
Multi-dimensional simulations will also allow a more de- 
tailed study of how the late-time relativistic outflow emerges 
from within the earlier non-relativistic wind and the differ- 
ent observational signatures associated with each phase of 
the outflow (a nalogous to studies of jets emerging from a 
host star; e.g.. iMorsonv et al.ll2007h . They also allow inves- 
tigation of the stability of the confining envelope, which is 
subject to possible Rayleigh- Taylor fragment ation, since it 
is accelerated by the light weight bubble fe.g.. lAronsll2003l ). 
an effect likely to be of substantial significance at higher en- 
ergy injection rates. Such instabilities could be important 
for understanding the short timcscale variability observed 
in GRBs, as could intrinsic time variability in the magnetar 
wind, driven by, e.g., reconnec tion in the equatorial current 
sheet (|Bucciantini et ail [2006'l. Variability might also origi- 
nate inside the bubble due to the dynamics of the interaction 
of the win d with the progenitor ( as is seen in PWNe simula- 
tions; e.g., Del Zanna et al ] |2004T ). Finally, in addition to the 
asymmetric expansion of the entire bubble discussed here, 
observations and simulations of PWNe also reveal a mod- 
erately relativistic axial "jet" within the neb ula itself (e.g., 
IWeisskopf et al.ll2000l ; iDel Zanna et al.ll2004h ; this may be 
dynamically important for proto-magnetar bubbles as well. 
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APPENDIX A: SELF-SIMILAR SOLUTIONS 

It is well known that the ID thin shell equations admit self- 
simi lar solutions, whe r e the time dependence can be factored 
out. lChevalier fc Lud (|l994T ) have shown that equations ([8]|- 
()13[) can be solved in the simplified case of a wind blowing 
inside a wind, when the shell expands at a constant speed. 
In this case it is possible to reduce equations (|8)l- (|13|) to a 
system of ordinary differential equations. Here we extend 
their treatment by considering the expansion of a hot mag- 
netized bubble inside either a stationary medium (as in the 
case of the ISM or the roughly static envelope of a progenitor 
star) or freely expanding ejecta (as in the case of PWNe). 
We derive the conditions under which self-similarity applies 
and discuss how the elongation of the bubble scales with 
its magnetic energy content. We limit our investigation to 
the case of an isotropic outer medium, even though self- 
similarity ho lds in the more general case of fixed latitud inal 
dependence i|Kahn fc Westll 19851 . ICampbell et al.ll2003l L 

If one assumes a cold unmagnetized outer medium with 
a power law density profile p oc r a , that p n oc t x and that 
Rpoiar/H is constant, where R po i a r is the radius of the shell 
along the axis, one can search for a solution of the form 

R = R{6)t 
v± = v±(9)t s 
wy = v n (9)t s 
a = a(8)t v 

Taking 5 = /3- 1, 7? = P + ap, and x = 2/3 + a/3 -2, the time 
dependence can be eliminated and the problem reduces to 
a system of ordinary differential equation in the 9 direction. 
Neglecting the outer pressure and magnetic field, and the 
inner ram pressure, and assuming a stationary outer medium 
(a similar derivation can be obtained for freely expanding 
ejecta), equations (|5))- (|13p can be written as 



idR 

Pa = p ((3R cos£) - 2/?o- - 
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(A3) 



a[v ll+ pRsmt]^ [£(„„) + PRco S £ (l + ||)] , (A4) 

where the variables are now only function of 8, and pi + 
Bl /8tt = pnf(Rs'm8/H). Note that these equations are not 
easy to write in non-dimensional form because, as we will 
discuss below, the characteristic length-scale is an eigenvalue 
of the problem. Only in t he case a = —2 : is th e solution 
relatively straightforward (| Chevalier fc Luolll994l ). 

To check the validity of the self-similar solution one 
must verify a posteriori that R po iar/H does not change in 
time and that the pressure p n scales as a power law. We 
first consider the scale height H . As discussed in |4] the 
evolution of the bubble can be described as a series of in- 
finitesimal adiabatic expansions followed by energy injection 
at fixed volume. It is easy to show that, in the relativistic 
case, when the adiabatic index is 4/3, entropy and magnetic 
flux conservation imply that the ratio R po iar/H is indeed 
constant during self-similar expansion. On the other hand, 
as long as the ratio of the Poynting flux to the total luminos- 
ity, E m ag/Etot, is equal to the ratio of the magnetic energy 
to total energy in the bubble E mag / Etot, the height scale H 
is constant, at fixed shape. 

Despite the fact that the thin-shell equations conserve 
only momentum, it possible to conserve energy as well by 
requiring the correct temporal scaling, in which case the 
total energy in the system increases only because of injection 
{Etot)- The total energy in the shell is: 



-cr(u 2 + v\)dS — E s ht 



5f3+cfi 



(A5) 



where 5* is the shell surface and E an is a constant that de- 
pends on the shell structure. If the ratio Rpoiar /H is con- 
stant then the total energy in the bubble is 



e(r, z)dv — p„Kt 



3/3 



(A6) 



where V is the volume of the bubble, K is a constant de- 
pending on the shell shape and the ratio R po iar/H, and e is 
the energy density in the bubble. If p„ oc t x then the total 
internal energy scales as the shell energy. Thus one can write 
the total energy in the system as 



E tot (t)= / Etotdt oc t sf3+a0 - 2 . 



(A7) 



This relation shows that as long as E to t is a power-law 
in time, self similar solutions are self-consistent. Equation 
(|A7|) together with the ratio E mag /Etot determines how 
the self-similar shape is related to the injection properties, 
namely the total luminosity and Poynting flux. Once a and 
Rpoiar /H are given one can search for a solution of equations 
(A1)-(A4). 



AO. 1 Numerical Methods 

In the case studied by IChevalier fc Luol l| 19941 ). the shell 
moves at a constant speed, P = 1, and the system of or- 
dinary differential equations can be easily solved as an ini- 
tial value problem, assuming regular conditions on the axis. 
Given that the equations are singular on the axis, it is nec- 
essary to expand them in Taylor series and solve for the 
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coefficients. In more general cases, the shell is subject to ac- 
celeration and /3 7^ 1. This changes the nature of the equa- 
tions and makes the problem much harder to solve. If one 
tries to expand the equations in a Taylor series on the axis, 
which works for the case j3 = 1, one finds that terms of order 
0(8 n ) depend on terms of order 0(8 n+1 ), and the series is 
open. This is typical of two-point boundary value problems 
(eigenvalue problems) where conditions at one boundary are 
not sufficient. However it is not obvious where and what kind 
of additional boundary conditions one should impose. One 
might enforce regularity on the equator, either by assum- 
ing a smoot h shape or zero p aralle l velocity. However the 
solutions of IChevalier fc Luol (|l994h show that discontinu- 
ities at the equator might arise, especially for larger value 
of Rpoiar/H. We have chosen to impose that the tangential 
transport velocity should be on the equator to avoid the 
formation of an equatorial ring where matter can accumu- 
late indefinitely, but this does not preclude other possible 
boundary conditions. 

We tried several methods for solving the steady state 
angular equations given by (A1)-(A4). These included stan- 
dard shooting techniques, using both explicit ODE solvers 
and an implicit integrator, standard relaxation methods, and 
a Fourier expansion of the angular equations. For a variety 
of reasons, we find that none of these methods was satisfac- 
tory for finding solutions, aside from the relatively spherical 
limit when R po iar/H <^ 1. Instead, we found that the most 
successful method was to use the time dependent equations 
(eqs. [8]- [13]) to evolve the system forward in time, imposing 
a fixed value for R po iar/H and a pressure p n oc t x , where 
the exponent x is the one expected for the self-similar so- 
lution. Self- similarity is then reached at late times after the 
initial transient dies away; we have verified that the shell at 
late times has converged to a self-similar shape, with only a 
few percent error in radius. It is, however, known that ac- 
celerated shells are subject to corrugational instabilities (the 
thin-shell instability). We have included an artificial viscos- 
ity to suppress this instability and focus on the overall evo- 
lution of the bubble. Instabilities may indeed be important 
for the evolution of the shell but we leave a study of them to 
future work. In each case, we have determined a posteriori 
that the chosen value of viscosity was small enough not to 
affect the shape of the bubble. This method for determin- 
ing the self-similar solution worked up to R po iar/H ~ 10. 
More elongated bubbles were difficult to investigate unless 
one started with initial conditions very close to the desired 
self-similar solution. Otherwise the deviation from the cor- 
rect self-similar solution was effectively a large-amplitude 
large-scale initial perturbation; in this case, we found that 
the corrugational instability could not be suppressed with- 
out increasing the viscosity to a point where it modified the 
overall shape of the bubble. 

AO. 2 Results 

In Figure lATI we plot the shape of the shell for three different 
values of R po iar/H, and for a = 0, — 1, —2 in the case of a 
constant injection luminosity (/3 = 3/ (5 + a)). For a = and 
— 1, the solutions are found by the methods described in the 
previous subsection and are r estricted to R m iq r /H < 10. 
The a = — 2 case is that of IChevalier fc Luol (l994) for 
which the angular equations (A1)-(A4) can be straightfor- 



wardly integrated. In this case, we compute solutions up to 
R po iar/H « 100, for comparison to the time-dependent so- 
lutions in the main text. 

Figure IA1I shows that the elongation of the bubble in- 
creases with increasing R po iar/H, and for smaller values of 
a. In the case of a constant density outer medium, the shell 
appears to be regular both at the pole and at the equa- 
tor, while in the case a = — 1 and —2 for large values of 
R po iar/H, a cusp is formed in the equatorial plane, where 
matter tends to accumulate. We have verified that all quan- 
tities scale according to self similarity. 

Unlike in the case studied bv lChevalier fc Lud l| 19941) , in 
which the surface density monotonically increases from the 
pole to the equator, for a = and small values of R po iar/H 
we find that the surface density reaches a maximum at in- 
termediate latitudes. 

Figure K2l shows a summary of the results of our simu- 
lations for the case of constant injection power into a static 
outer medium (typical for wind bubbles): shown are the ratio 
of the magnetic energy to the total energy and the elonga- 
tion of the bubble, as functions of R po i ar /H. The elonga- 
tion scales approximately as -J R po iar/H and, as expected, 
is larger for steeper density profiles. The dependence on the 
density profile of the outer medium is, however, quite weak, 
with less than a factor of 2 difference between a — and 
a = —2. For the magnetar problem considered in this paper, 
this suggests that the outer density structure of the progen- 
itor star should not significantly effect the collimation of the 
bubble (this is consistent with the fact that the asymmetry 
of the bubble was only a weak function of the progenitor star, 
as discussed in . Figure IA2I also shows that in all cases 
the ratio E ma g/Etot inside the bubble is relatively small, 
and tends to be smaller for smaller values of a. This demon- 
strates that strong elongation does not require magnetically 
dominated bubbles, which is consistent with our result for 
the evolution of magnetized bubbles inside GRB progenitors 
in g] 

We also carried out simulations for a bubble inside 
freely-expanding ejecta with velocity v oc R/t and density 
p oc r a t~ a ~ 3 . We find that for a = —2, the time- dependent 
solution never reaches a self-similar solution, but instead the 
evolution always depends on the initial conditions; for this 
reason, we do not show solutions for a — —2. This because 
in the case a — — 2 if one assumes the expected power law 
temporal evolution for the internal pressure, it is easy to 
show that the ratio between the inner pressure and the ram 
pressure of the outer medium is constant, and does not de- 
pend on the radial evolution. In this sense, the bubble never 
relaxes to self-similarity. The case of freely-expanding ejecta 
is relevant for a central source which is not energetic enough 
to produce a bubble which effects the dynamics of the stellar 
explosion on short time-scales, but instead creates an ener- 
getic PWN inside the expanding SN ejecta. This could be 
relevant to young PWN generally and also to lower energy 
transients such as X-ray flashes (see §5). The results for the 
elongation of the bubble and its dependence on the magne- 
tization are shown in Figure DOl they are reasonably similar 
to the results shown in Figure POl 
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Figure Al. Self-similar shapes for magnetically driven bubbles expanding into a static ambient medium with p oc r a . Results are shown 
for three different values of R po i ar / 'H . Results for a = —2 are based on the model of lChevalier fc Luol lll994Tl while for a = and — 1, 
they are based on time dependent thin-shell simulations that reach self-similarity at late times. 
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Figure A2. Properties of magnetized bubbles expanding into a 
stationary medium with a power law density profile with p oc r a , 
for a = (squares), a = —1 (triangles), and a = —2 (diamonds). 
The a = — 2 results are based on solving the steady state angu- 
lar equations (Af)-(A4), which can be solved for larger values of 
Rpolarl H than the a = and —1 cases, which are based on solv- 
ing the time dependent thin shell equations. E mag / Etot is the 
ratio between the magnetic and total energy in the bubble. The 
elongation is defined as the ratio between the length of the bubble 
along the axis and the maximum extent of the bubble away from 
the axis. 
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Figure A3. Same as Figure lA2l but for the case of expansion into 
the freely expanding ejecta of a supernova remnant: v oc R/t, 
p oc T a t~ a ~ i ; results for a = (squares) and a = —1 (triangles) 
are shown. 
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